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We introduce the idea of shape parameters to describe the shape of the pencil of rays connecting 
an observer with a source lying on his past lightcone. On the basis of these shape parameters, we 
discuss a setting of image distortion in a generic (exact) spacetime, in the form of three distortion 
parameters. The fundamental tool in our discussion is the use of geodesic deviation fields along a null 
geodesic to study how source shapes are propagated and distorted on the path to an observer. We 
illustrate this non-perturbative treatment of image distortion in the case of lensing by a Schwarzschild 
black hole. We conclude by showing that there is a non-perturbative generalization of the use of 
Fermat's principle in lensing in the thin- lens approximation. 

I. INTRODUCTION 

Before we describe the ideas and content of this work, we wish to give a brief discussion of certain ideas concerning 
lensing that are really quite obvious but are often overlooked and or confused. It is absolutely clear that any source 
of light emits many rays that travel in many directions and get distorted by all the different available lenses. But for 
an observer, of all those different rays, the only ones of relevance are the rays that enter the observers "eyes;" all the 
others are "lost" to that observer and play no role at all in his or her treatment of lensing. The conclusion from this is 
that at any one time, it is the past lightcone of the observer that is the object to be studied and understood. It is the 
bending and distortions of the rays and pencils of rays lying on this past cone that constitute the basic elements that 
make up lensing theory. When this past cone intersects a source, that source is then seen by the observer; if it is a 
point source then ideally only one ray reaches the observer, if it is a small but finite source then a pencil of rays, lying 
on the past cone, reaches the observer. In this later case the pencil can be described by the geodesic deviation vectors 
(or Jacobi field) associated with just one ray. The shape of the pencil can and will, in general, undergo changes due to 
the presence of gravitational field along the path of the one ray. It is this change in the shape of the pencil, between 
the source and the observer, that constitutes image distortion. This is the point of view that we adopt throughout 
this paper - though the details are very different, the point of view concerning the essential role of the past cone 
applies equally well to the thin lens approximation (the standard approach) and to the non-perturbative approach to 
lensing. 

It is often the case in gravitational lensing that the distortion of extended images provides as much information 
to the practicing astrophysicist as the location of the images themselves |^ . In previous work [||-|| , the theory of 
gravitational lensing has been developed and studied from a first principles point of view, where lensing is considered 
non-perturbatively with no restrictions on the strength of the gravitational field. This work is non-perturbative in the 
sense that the key equations of lensing theory, the lens and time of arrival equations, are derived by integrating (as far 
as possible) the null geodesic equations of an (in principle) exact metric. This point of view is in contrast to the more 
standard work in lensing based on the thin lens approximation where lenses "live" in some fixed background spacetime. 
In the present paper and its companion Q , we turn our attention to the non-perturbative study of distortion of images 
due to lensing. 

Images are distorted by a lens in multiple ways. First, an image may be magnified or demagnified, relative to 
what the image would look like "in the absence of the lens." Magnification is due to the focusing of the lightcone 
by the spacetime curvature, and the study of magnification is closely related to catastrophe optics and singularity 
theory []6|,[7|jlC|-]r^ . The magnification of images has proved to be an extremely useful property of lensing systems 
because it often allows for the observation of very distant sources which may not have been observed without the 
magnification. 

Perhaps more importantly, an extended image may be stretched along some axis, rotated, or made generally larger 
or smaller relative to an "unlensed image." Colloquially, we say that the stretching along an axis and rotation of an 
image is related to the "shearing" of the pencil of light and that an overall change in size is due to the "convergence" 
(or alternatively, "expansion" or "divergence" ) of the pencil. These effects are very important in lens modeling, as 
the distortion of an extended image into an arc or of identifiable parts of an image (such as the radio jets of a QSO) 
provides very useful observational constraints on the lens model. 
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In general relativity, the shear and the convergence are precisely defined quantities which have values at each 
point along a geodesic. These quantities, often refered to as the optical scalars, have been widely studied in the 
relativity community. Hence, the general use of these terms in lensing, where an image - observed here - of a source 
- over there - is said to be "sheared" can be a little misleading. 

One goal of this paper and its companion is to propose a unified framework for the non-perturbative discussion of 
the distortion of images which can be applied in strong-field regimes. In this paper, we introduce three quantities, 
the shape parameters, describing the generic appearance of an elliptical image or source and show how the shape 
parameters are computed from connecting vectors. We also introduce three distortion parameters which correspond 
to changes in the shape parameters between the source and observer. In the following paper, we relate the shape and 
distortion parameters to the optical scalars. 

This paper begins by reviewing the standard (thin-lens) approach to the distortion of images as presented in psl. 
In Sec. II, we introduce the key idea of using the inverse Jacobian of the lens equations to map source curves (the 
outline or boundary of sources) into image curves. We define the shape parameters at the end of this section and 
briefly mention the quantities which are referred to as the "shear" and "convergence" in the thin lens approximation. 
(The full discussion of these quantities is left to the second paper, which is devoted to the relationship between the 
shape parameters and optical scalars.) In the third section of this paper, we present the non-perturbative version of 
the shape parameters and define the distortion of images. This section relies heavily on the use of connecting vectors 
or Jacobi fields (we follow the lead of other authors in the use of Jacobi fields in lensing, in particular |8||]). The 
distortion of images in Schwarzschild spacetime is discussed at the conclusion of this section. Finally, we provide a 
different perspective on the subject in Sec. IV which wc feel clarifies several issues regarding the role of the choice of 
coordinates. 

An important caveat must be stated. In a generic lensing situation, if the past lightcone is continued arbitrarily 
far into the past, caustics will always be found. We, in this work, are making the assumption that our source that is 
being imaged does not intersect a caustic curve. The morphology of image construction for sources which lie across a 
caustic is well understood but relatively complicated. The issue of distortion of images in this case is something for 
further study. 



II. IMAGE DISTORTION VIA THE THIN LENS APPROXIMATION 

In this section we summarize the approach to image distortion in standard gravitational lensing according to p5| . 
The assumptions of the thin lens approximation are that the source, lens, and observer are at rest in the background 
spacetime, that the gravitational fields are weak, that all angles are small, and that the mass distribution of the entire 
lens can be collapsed into a single two-dimensional "lens plane" and is described by a surface mass-density, k(^). As 
is shown in Fig. |l|, light rays are assumed to travel along the null geodesies of the background spacetime from a point 
in the source plane at S" to a point in the lens plane (the plane containing the lens and perpendicular to the optical 
axis) at /. At /, the ray is instantaneously bent by a bending angle, a, and travels to the observer O along another 
background geodesic. A goal of the thin lens approximation is to find a set of equations, called lens equations, which 
represent the path of the light ray as a mapping between the lens plane and the source plane. 

The lens plane is described by the two dimensional vector ^, and the source plane is described by the two dimensional 
vector ff. It is convenient to introduce two scaled vectors 

C V 
X — — , and y = , (1) 

which describe the lens {x) and source {y) planes. Here Ds and Di are the angular-diameter distances between the 
observer and source and observer and lens. 

The dynamics of the thin lens approximation is introduced through a deflection potential 

Mx)^- [ dx'^K(x')\n\x-x'\, (2) 

where k(x) is the surface mass-density, satisfying the two-dimensional Laplace equation 

VV(^) = 2k(x). (3) 
The lens mapping from the lens plane to the source plane, or from x to y, is a gradient mapping: 



2 



y {^x-x-^{x)] , (4) 



y ^ X — 'S/il){x). (5) 

Equation is the lens equation of the thin- lens approximation under the scaling defined by Eqs. (|l|). 

An alternative perspective is to derive the lens equation, Eq. , by applying a stationarity condition to a Fermat 
potential. If we define 

^{y,x) ^'^{x-yf ~tl){x), (6) 
as a Fermat potential, the gradient mapping is equivalent to the condition 

V,0(y, x) = 0, (7) 

i.e., the lens equation, Eq. (||), follows from Format's principle. (We generalize the idea of a Fermat potential in Sec. IV 
to a general lensing scenario, noting that other authors have also studied this problem extensively fl^.) Note that 
the second derivativcof (j){y, x) with respect to x is equal tothc derivative of y with respect to x as defined inEq. (|). 
The distortion of images can be studied via the Jacobian of the lens mapping 

A=^ = (8) 
dx dxdx 

We see that the Jacobian matrix is symmetric. This symmetry exists, essentially, only in the special coordinates used. 
We will see in Sec. IV how this result arises non-perturbatively. 

Using the lens equation, Eq. (||), we see that the Jacobian matrix has the form 

where subscripts denote derivatives with respect to the components of a; = By defining 

7l(^) = ^(V^ll - "022), 

72(f) = -012 = V'21 (10) 

and using Eq. (|^), we can rewrite the Jacobian matrix as 

A=f l-'^-^i -^l V (11) 

V -72 1 - K + 71 y ' 

where we have dropped the explicit x dependence in k{x), 71(2?) and 72(2;). 

It has become customary in the lensing community |15 1^ to refer to k as the convergence, and 7 = ^77^ + 71 as 



the shear . In this and the subsequent paper ||], however, we will reserve the names of convergence and shear for the 
optical scalars of a null geodesic congruence. We show in our subsequent paper that (modulo terms involving Ds and 
Di) K and 7 are essentially the values of the convergence and shear of the observer's past lightcone at the lens plane 
on the source's side. 

In terms of k and 7, the invariants of A are 

detA={l-Kf~-/^ (12) 
tr A = 2(1 - k) (13) 
a± = 1 - K ± 7 (14) 

where a± are the eigenvalues of the Jacobian matrix. When det A = 0, the lens mapping is singular and the source 
lies at a caustic of the past cone. The quantity /i = (det A)"^, referred to as the magnification, is the factor by which 
the image's brightness is magnified with respect to the unlensed source. 
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Another important use of the Jacobian is to help understand the changes in the shape of the image of a small 
source. To see this, we note that a (suitably scaled) connecting vector, X, at if in the lens plane, is mapped into a 
(suitably scaled) connecting vector, Y , at y in the source plane, by A, while maps source plane vectors back to 
the lens plane: 

Y = AX. (15) 

For instance, if a small circle of radius i? at ?/ in the source plane is described by 

Y{t)=y^Rs{t) (16) 

where s[t) = (cosi, sint), an image of this circle is an ellipse, in the lens or image plane, described by 

X{t) = x + RA-^s{t) (17) 

Hence, the Jacobian matrix. A, and its inverse, A^^ , which is referred to as the magnification matrix, play a large 
role in the description of the images of small sources. 

The semiaxes of the ellipse described by Eq. ( p^ ) must be parallel to the principal axes of the magnification matrix 
A^^, which are the same as the principal axes of A. This is because there are two vectors s'{t±) = {cost±,smt±), 
which are parallel to the two eigenvectors of A and are mapped into themselves by A and A~^. These two directions 
undergo the maximum stretching and squeezing because the eigendirections of a symmetric linear map are those which 
have the greatest and least stretching. Hence, the lengths of the major and minor semiaxes are 

L± = TT^—r (18) 
|1 - K=F7l 

The locations of these axes are given by the angles, t±, which satisfy 

(19) 

72 y \72/ 
i.e., in the directions 

X{t±) = x + RA-^s{t±). (20) 

This result is obtained [via extremizing {X{t) — xY as a function of t] by finding the angles, t, for vectors of the form 
of sit) = (cos i, sin t) which are the eigenvectors. 

We are now in the position to define our shape parameters. The first shape parameter is the ratio of the lengths of 
the semimajor and semiminor axes, TZ. This ratio is given by 

\l-K-^\ ^ ' 

in the thin-lens approximation. The second shape parameter is the area of the ellipse factored by tt, or the product 
of the two semiaxes, which we refer to as the area of the ellipse (since the area is proportional to the product of the 
semiaxes with a factor of tt). For the thin lens approximation, the (factored) area is given by 

AeeL+L_ = - - = R^det{A-') (22) 

(1 — ny — 7^ 

The final shape parameter is the orientation of the major semiaxis, 5. By convention, we define 5 to be (5 = t+. 

Thus, there are three quantities of interest in describing the image of a circular source: the (suitably scaled) area 
of the image A, which encodes the overall scaling of the image compared to the source; the ratio of the semiaxes TZ 
which is a measure of the distortion up to scale; and the orientation of the semimajor axis 5 = t-^-. This definition of 
the orientation does not have a definite meaning without a reference axis or when referred to a circular source but it 
can be assigned a meaning for an elliptical source when the discussion is generalized slightly as will be done later in 
this article. We refer to these three quantities as the shape parameters. They will be defined non-perturbatively in 
the next section. 
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III. NON-PERTURBATIVE IMAGE DISTORTION 



We are interested in understanding image distortion within general relativity with no approximations. We refer to 
this approach to image distortion as a non-perturbative approach, to distinguish it from the standard approach which 
is based on linearized perturbations around flat space or an isotropic cosmological model. 

To this effect, we need two things. First, we need a non-perturbative counterpart of the lens mapping. Once the 
lens mapping is given, the Jacobian matrix, or its inverse, the magnification matrix, can be obtained. Second, we 
need non-perturbative definitions for image and source shapes, and distortion parameters, which encode the changes 
between the shapes of the source and the image, in terms of the magnification matrix. 

A. Lens mapping 

At any fixed time the observer "sees" each source (on his or her past lightcone) projected onto the observer's celestial 
sphere. By introducing spherical coordinates, {9, 4>), for each ray, via this projection, the observer encodes a mapping 
between an image location, {0,4i), and a source location. A formal lens mapping can be obtained by the following 
considerations, summarized from [^,^. Similar considerations have been studied recently by other authors ||l^ , pj| . In 
particular, see for exact lensing in Kerr spacetime. 

In a four dimensional Lorentzian space-time, {M, gab{z"')), with local coordinates we assume that an expression 
for the past lightcone of the observer, obtained by integrating the null geodesic equations, is available. It has the form 
of four functions of three parameters of the light cone at an observer's fixed time, r: 

z-^F-{z',{T),s,e,cl,). (23) 

Here z° represents a point on the past lightcone of the observer's worldline, Zq{t) at time r, reached at an affine 
length s along the null geodesic with direction (6*, The affine length s has the value zero at the observer and is 
scaled so that gabF°'ZQ = 1, using the notation ' = Thefunctions can be obtained by integrating the geodesic 
equations, F" + Tl^P'^F'^ = with the null condition 

gahF'^F'' = 0. (24) 

The coordinates are arbitrary, but they can be chosen so that one of them, say z°, is timelike and the remaining 
three are spacelike, i=l,2,3. In this case, for a ~ 0, Eq. (^3|) represents the time of departure of a light signal from 
the source to the observer. On the other hand, Eq. ( p3| ) for a = 1,2,3 is equivalent to a lens mapping between the 
angular position, (z^,z^), say, of a source of light and its image angles (0,0) on the observer's celestial sphere, plus 
a prescription of a radial coordinate, z^, which, in principle, can be obtained from independent observations. Ideally, 
if the radial coordinate is known as a function of the affine parameter s, then, in principle, the parameter s can be 
eliminated from F^{s,6,(j)), for A = 2,3, in terms of {z^,9,(f>), in which case a precise analog of the standard lens 
mapping is obtained of the form = {z^ , , (j)) . For our purposes, it is not essential to eliminate s, which may 
not even be possible in general (see [Q for the feasibility of such an elimination in the simple case of a Schwarzschild 
spacetime). The non-perturbative lens mapping z"^ = G^{z^,6,(j)) has been worked out explicitly in the case of a 
Schwarzschild spacetime Q and compared to standard lensing. 

Aside. Our non-perturbative version of a lens mapping differs from the standard lens mapping, in addition to being 
non-perturbative, in the sense that we do not use explicitly the Einstein equations - we simply consider any Lorentzian 
metric gab{z°') and its null geodesies to construct the observer's past cones and the resulting lens equation. In most 
thin-lens treatments of lensing the linearized Einstein equations (in the form 'S/'^ip{x) — 2k{x)) are used explicitly 
from the start in the construction of the lens equation. 

B. Jacobian of the Lens Mapping 

Consider the situation where there is an extended source at some distance from an observer. The source's worldtube 
intersects the observer's past lightcone, and we assume that the intersection occurs in a region of the past lightcone 
where there are no caustics. See Fig. |^. The intersection determines the source's visible shape. The pencil of rays 
that joins the shape of the source to the observer "carries" the shape of the source into the shape of the image, on the 
observer's celestial sphere. For the lensing purposes, thus, an extended source consists of a collection of individual 
point sources, all of which are connected to the observer by individual null geodesies. In general, the image of an 
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extended source must be obtained by imaging each individual point via the non-perturbative lens equation. However, 
if the sources are "small," we can, in the spirit of the previous Section, describe sources by connecting vectors along 
one of the null geodesies (a "central" one), with image directions {9, (j)), that joins the source to the observer. Jacobi 
fields (or connecting vectors) are solutions to the geodesic deviation equation along that fixed geodesic. In the case 
of null geodesies, the space of Jacobi fields is two-dimensional and can be obtained in our approach directly from our 
non-perturbative lens mapping Eq. (^3|) in the following manner. 

We note that there are three natural vectors that can be associated with either the past lightcone or equivalently 
with the lens equation, Eq. (|2^). The first one is the tangent to the geodesies in the lightcone: 

r EE ^ = (25) 

OS 

The remaining two are the Jacobi fields 

. rj ' f)z^ 

Mf^— , M^^{sm0)-' — . (26) 

These two vectors are, in general, linearly independent and connect neighboring geodesies at a fixed value of s, so 
they are two linearly independent Jacobi fields (except at a caustic). They are orthogonal to the tangent vector 
because they lie on the lightcone: gabM^fi — gabM2i'' — 0, and they are spacelike. Though any other Jacobi field is 
a linear combination of (M°,M2) with constant parameters a, (3, 

Z" = aM^ + f3M^, (27) 

we refer to (Mf ,Mf) as the "natural" Jacobi fields. 

Two other important vector fields (not Jacobi fields) associated with any null geodesic are the pair of spacelike, 
orthonormal, parallel propagated vectors that are normal to They provide a means to refer to directions that are 
"the same" at the orthogonal tangent spaces at the source and at the observer. Analytically they are given by 

ff-Vbe? = £''Wbe^ = 0, (28) 
ei.£ = 62-^ = 0, (29) 
ei- ei = 62- 62 = -1, (30) 

61-62=0, (31) 

where v-w denotes the scalar product between two vectors, gabV°^w^ ■ 

The natural connecting vectors, (Mf , ), in general, will have transverse components along the space spanned by 
(ef , 62), and a longitudinal component along the geodesic congruence t"'. The component along has no significance 
since it corresponds to reparametrizations of the null geodesic. Jacobi fields thus have two essential transverse 
components. Two linearly independent Jacobi fields can thus be thought as two two-component vectors that are 
obtained by derivation from the lens mapping. Explicitly, this can be done by the projection of {Ad^.M!^) into 
{el, 62), thereby defining the matrix J(s) 

{Ml, M-) ^ [j{s)] EE j\{s) = M!^el (32) 

At the observer (i.e., at the apex of the cone), the connecting vectors (M",]^^) vanish, so that J(0) — 0, and, 
furthermore, in the limit as s — > we have that 

lim -J{s) = I (33) 

s^O S 

where I is the identity matrix. This last result is quite important and comes from the fact that although the connecting 
vectors vanish at s = 0, their angular size has been chosen to be unity. (At each value of s close to the observer, each 
connecting vector subtends an arc from the observer's viewpoint. The connecting vector's arc is proportional to s, 
therefore the angular size of the connecting vector at the observer is well defined, and equal to the derivative of the 
connecting vector with respect to s.) 

We adopt the following notation. When the matrix J{s) is evaluated at the source position, denoted by s = s*, it 
is called the Jacobian of the lens mapping and is denoted by *J. In addition, we denote any quantity that is defined 
or evaluated at the source by placing a * in front of it, e.g., {*M^, *M^) = {M^{s*),M^{s*)). 
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Notice that, in a weak field-thin lens regime, our Jacobian *J is related to the Jacobian A of the thin-lens theory 
via 

*T 

A^- (34) 

We now put a source at s* with an elliptical shape and study how it is mapped back to the image space, the 

observer's celestial sphere. Though it is possible to do this construction with (Mf (s), Aff (s)) (or,equivalently, with 

— j 

J,;(s)) it turns out to be considerably simpler to use a pair of Jacobi fields, (Mf,M|), constructed from a hnear 
combination of the Mf that are unit orthogonal vectors at s = s* , i.e., so that Jj{s) = Mfe^ is the identity matrix 
at s = s*: 

Ji{s*)^*Mrei=Si. (35) 
It is easy to see that the linear combination must have the form 

Mf = i*r')lM^, (36) 

so that 

M-e^ = (y~')^M-e^ ^ J^s) = (tT^)^ j'(.), (37) 

4(s*) = (tr')^tf* = 5f. (38) 

This is equivalent to assuming that we have chosen two connecting vectors and so that 

*Aff = *e1, (39) 
= *e^. (40) 



and hence 



In general, the two pairs (M{',Af2) and (6^,62), coincide at the point s* , but not at any other point along the null 
geodesic. See Fig. 0. 



For our construction in this article, it is a fundamental observation that, from Eqs. (pT\) and (33), as s — s- we have 



-J{s)^U \ (41) 
s 

i.e., in the limit as the null geodesic reaches the observer, J acts essentially as the inverse of the Jacobian of the lens 
map, namely, as our magnification matrix. We make extensive use of this result throughout this and the companion 
paper [§. 

Preliminarily to defining the shape and distortion parameters we point out that even though the pairs (e°, 62) and 
(Af^Aff) possess different meanings from each other, they can nevertheless be expressed as linear combinations of 
each other. Assuming a relationship of the form 

= a cos(A + v)M1 + (3 cos XM!^, (42) 
= a sin(A + v)M1 + (3 sin \M^, (43) 

and using the properties of e° and its dual basis, namely 

e^ei = 5l e^4 = 0, i'V.el ^ Q (44) 

after a straightforward calculation the four parameters (a, /?, v, A, ) can be obtained in terms of the scalar products 
Ml- Mi,M2 - M2 and Mi ■ M2 as follows: 



-A/2- 


M2 


C 


-Ml- 


Ml 


C 




Ml 


■M2 



(45) 
(46) 

cos ly = .,.1.,.., ^ 
VAfi-A'/iA/2-A'/2 
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where we use, for short, 

C = Ml ■ M1M2 • Af2 - (Afi • A^2)^ (48) 

and 

A^^^logf^V (49) 



2x/C ds ^\Nh-M2 

which fixes A up to a constant. The quantity C represents the area spanned by {M^^M!^)^ i.e., C — g'^'^UaUd with 
Ua = eahcM^Ml 

Notice that the values of (a, /?, A) at s* are determined to be 

We thus have the orthonormal paraUel propagated basis written in terms of the Jacobi fields. 



C. Shape and distortion parameters 



We now discuss the source parameters that enter into the considerations of image distortion. Specifically we consider 
small elliptical sources, which can be completely described by three parameters: the lengths of the two semiaxes, and 
the orientation of, say, the major semiaxis with respect to a reference direction. The two semiaxes can be equivalently 
encoded into the ratio of major to minor semiaxis, and the product of the two. The product of the two semiaxes is 
the area of the ellipse divided by the number tt. In order to avoid carrying the factor tt around, in this paper, we 
simply refer to the product of the two semiaxes as the area of the ellipse and denote it by *A. We denote the ratio of 
the semiaxes by *TZ, and the orientation of the major semiaxis with respect to the parallel propagated direction by 
*5. These are, by definition, our shape parameters evaluated at the source, s* . 

At the source location, we have an orthonormal basis, (*e'f , *e2), in terms of which an ellipse can be expressed as 

V(t) = Y^{t)*el + y^(i)*e2 = Y^{t)*M1 + Y^{t)*M^ (51) 
where the components Y^{t), t being the curve parameter, can be given by 



\ _ ( cos*(5 - sin*(5 \ / VWR 

V2 — n 



cosi 



/ZTm? { COS i cos*(5 - sin t sin*^ \ , . 



Y^ J y sm*S cos*(5 J \ ^/*A/^ J \ sint 

IZ cos t cos*(5 — sin t sin*S 
IZ cos t sm*S + sin t cos*S 

As indicated, this arises from first taking a parametrized circle (the vector (cost, sini)), then distorting into an ellipse 
of major semiaxis V i4*7?. (the product of the semiaxes times the ratio of major to minor is the square of the major 
axis) and minor semiaxis y/*A/*TZ (the product of the semiaxes over the ratio of major to minor is the square of the 
minor axis) and finally rotating through an arbitrary angle *S. 

Remark: One should not confuse the curve parameter t with the polar anglecj) of each corresponding point on the 
curve relative to the direction *e". The polar angle (p is related to t at the source by tan((p ~ S^,) — tant/*R. 

In order to obtain the image of this source, we "follow" each point of the ellipse back towards the observer using 
the Jacobi fields Mf(s), i.e., via 



Z''{s,t) = Y\t)M^{s) + Y^{t)MS{s) 
The vector field Z'^ is a Jacobi field for every value of t, therefore it describes the image of the source near s ^ 0, the 



■ (i*n cos t cos*(5 - sin t am*5)Aq + {*R cos t s\tl*5 + sin t cos*5)M^^ . (53) 



observer's location. On the other hand, it agrees with our expression for the source at s* , because of Eqs. (pS9|)-(40). 
Projecting the components of Z"'{s, t) in the ef directions we have 

Z\s,t) Z^{s,t)el Y^it)Mns)e: = J){s)Y'^{t), (54) 
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the matrix J thus pulls the source's shape back towards the observer. We then have the projected image of the source 
on the observer's celestial sphere given by 

X\t) = hm -J]{s)Y^it) = i*r')]Y^t) (55) 

or 

X{t) = tT V(i). (56) 

This is the non-perturbative equivalent of Eq. ( [T7| ) of Section II. Notice, however, that we have not scaled the 
connecting vector at the source. (The connecting vector at the observer must be scaled by s in order for it not to 
vanish.) 

The image X'^{t) of the ellipse Y"-{t) is an ellipse as well because a non-singular continous linear transformation of 
an ellipse is also an ellipse. In order to read off the image's elliptic parameters, we start by extremizing the distance 
Z ■ Z as a function of the parameter t. Directly from Eq. ( [S^ ) we have 

*A 

Z- Z = — {a + bsin t + csintcost) , (57) 

with 

a{s) = *7^2(^(cos*(5)2Ml•Ml + {sm*5)^ M2 ■ A'h + 2 sin*(5cos*(5A/i -Ma) , (58) 
b{s) ^ i(l -*;^2)(Ml•Ml + Ma-Ma) + ^(1 +*TZ^) cos{2*6){M2- M2 - Mi-Mi) - (1 +*7^2) sin(2*5)Mi -Ma, (59) 
c(s) = *7^(^sin(2*(5)(Ma•Ma - Mi-Mi) + 2 cos(2*(5)Mi -Ma) . (60) 
Extremization {d{Z-Z)/dt — 0) yields 

tan(2i) = -7. (61) 


There are, thus, four values of t in the interval [0, 27r] that extremize the distance to the center, and they differ by 
7r/2. Using Eq. (|6l|), we can obtain, for Z-Z, 

Z-Z=^(2a + b±Vb^ + c^). (62) 

The lengths of the major and minor semiaxes, and L_, respectively, are thus 

1/2 



fJA_ ' ^ 



^-=\jk (-2«-&-V62 + c2)) ■ (64) 



The ratio of the semiaxes is thus 



1/2 



7e4^" + ^-^) (65) 

where a,b,c are explicit functions of the scalar products between the two connecting vectors Aff and M2 given by 
(§)-(©• 

The product of the semiaxes is what we refer to as the area of the ellipse 

A=^{M^ + b)-c'f'. (66) 



With the parameters a, b and c given by Eq. (pSD, Eq. ( p9D and Eq. ( pOj ) this becomes 

1/2 



A = a(^Mi ■ M1M2 ■ Ma - (Ml • Ma)^) (67) 
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as is expected. 

In order to find the orientation 5 of the ellipse with respect to we need to express in the basis (e°, 62). This 
is done by inverting Eqs.p^)-(|4^), obtaining 

, cos AeS ~ sin Ae? 
a sm v 

j^ja ^ sin(A + t/)e^ - cos(A + iy)e^ 
^ P sin ' 

The orientation S is then given by tan (5 = (Z- ei)l(Z- ei), or 

^ ^ acos(A + i')(7?.* costsin(5* + sintcosJ,) — /3cosA(7?., costcosJ, — sinisinJ*) 
/3sinA(7?., cosicosJ, — sinisin(5,) — asin(A + v){Jl^ costsin(5, + sintcos(5») 

For generic values of t this yields the polar angle that the points on the ellipse make with the direction e\. The 
orientation of the semiaxes is obtained by using, for the values i± of Eq. (^l]). It is straightforward to check that 
the values of t± given by Eq. (|6l|), which differ by 7r/2, yield values of 8 that also differ by 7r/2, as expected for an 
ellipse. [Other authors ||2^] have studied the rotation of the image of an elliptical source but in a different physical 
context (the appearance of the polarization of radio emission by sources at high redshift).] 

We have thus obtained our shape parameters 7?., A and b in terms of Jacobi fields, in Eq. (|65|), Eq. (|67| ) and Eq. (|7C|). 

In the limit as s — > we can obtain, from these shape parameters, the values of the elliptic parameters of the 
image. Notice that both and (5 arc well defined and finite as s 0, because they are ratios of connecting vectors 
that vanish at the same rate (like s). Therefore, if we denote the ratio of semiaxes of the image by 7?./, and their 
orientation by (5/, their values are simply the values of 7?. and i5 at s = 0: 

7^/ = 7^(o), (71) 

81 ^m- (72) 

The area of the image, however, is not a well-defined concept, since the image lies on the celestial sphere of the 
observer, which has no definite radius. Furthermore, the shape parameter that gives the area of the pencil or rays, ,4, 
necessarily vanishes at the observer's location (the apex of the lightcone). Nevertheless, the solid angle of the image 
(i. e., its "angular" area) is a meaningful concept. In order to obtain the image's solid angle, we define the solid angle 
51(s) of the corresponding the pencil of rays: 

^{s) ^ ^ (73) 

from which the image's solid angle fJ/ can be obtained by passing to the limit as s — > 0, namely: 

f7/ = limf7(s). (74) 

s— 

We define now the distortion of the image with respect to the source as follows. We consider the image to be 
distorted if its semiaxes and orientation are different from those of the source. Thus (L_|_/L'j_), (L_/_L1), and ((5 — (5*) 
are indicators of distortion. The image will be said to be distorted if any of the following distortion parameters is 
nonvanishing 



- 1, (75) 



2 



2)5 = (5(0) ^*6 = 5i - *6. (77) 

These distortion parameters give a measure of "total" distortion through the gravitational field between a distant 
source and the observer. The distortion accumulates, however, infinitesimally along the null path. It is the object 
of our subsequent paper to find the relationship between the "infinitesimal distortion" through an infinitesimal 
displacement in the gravitational field, and the optical scalars (convergence and shear) of the lightcone. 

We conclude this section with an illustration of our approach to distortion in the case of a Schwarzschild spacetime. 
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D. Example: The Schwarzschild Case 



As an illustration, we consider now the case of gravitational lensing by a Schwarzschild black hole of mass m. The 
line element is 



ds' =\^1-^^ dt^ _ (^1 _ ^ j - r\d9^ + sinHd(t>'). (78) 

In this subsection only, we are using {6, (f)) as the regular spherical Schwarzschild coordinates, instead of as coordinates 
on the observer's celestial sphere. The null geodesies can be given U,|2^ in terms of two initial angles at the observer, 
denoted -0 and 7, representing the observer's celestial sphere, and an inverse radial coordinate / = {^/2r)~^ , in place 
of an affine parameter s: 



t = ta + 2 



dV 



- 2^/2mlp) - P{1 - 2V2ml) V2l'^{l - 2^ml') 



(79) 



I = /, (80) 

cos 9 = — cos 9q cos Q + sin 6*0 sin 8 cos 7, (81) 

^ ^ sin 00 sin 6*0 — tan 8 (cos 00 sin 7 — sin 00 cos 7 cos 6*0) ^^^-j 
cos 00 sin 9o + tan 8 (sin 0o sin 7 + cos 0o cos 7 cos 9o) ' 

where (^o, ^o, ^0, 0o) are the coordinates of the observer's location, Ip is the inverse radial distance of closest approach, 
defined as the smallest of the positive roots of 

sinV^p(l - 2V2mlp) - Zg(l - 2V2mlo) = 0, (83) 

and where 8(/, Iq, ^) is the following function: 

f'p dl dV 

8(Mo,^) = vr-2 / , - / , (84) 

•^'o ^ll{l-2^/2mlp) -r\l-2^/2ml) -^i ^ l^il - 2V2mlp) ~ l'^ {1 - 2V2ml') 
for positive values of 0. For negative values of we have 

8(?,?o,V') = ~e(/,/o,-0), (85) 

namely, 8 is an odd function of ■0. Equations (|79|)-(|8^) describe points on the past lightcone of the observer which 
lie at an affine distance s larger than the afhne distance of closest approach, so they can be interpreted as source 
points behind the lens. Eqs. ([79|)-(p^ are the equivalent of Eq. (p3|), and Eqs. (pl|)-(p2|) represent strictly the lens 
mapping between the image's angular location ("0,7) and the source's angular location (0,0). By derivation of the 
lens mapping we define the starting connecting vectors 

^ ~ \di^' dip' dij' di^ ) \dip' 'di^'di^)' ^ ' 

^7^.(|i,«f.|^U(o,o, ,87, 



97 97 dj J \ 97 97 



where t, 9, are functions of (/, ip, 7) given by Eqs.(|79[), (pl|), (p2[), and the partial derivatives are taken at fixed value 
of I. See Figure ^. (Mf , ) are actually different from the "natural" connecting vectors introduced in Section III in 
that (MfjMf) are not normalized in a manner that would ensure their angular size to be 1. We have 
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de 



do 



sin Q sin 6^ sin 7 

s6q — sin 6 sin 6*0 cos 7)^ ' 
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^7 ^1 — (cos 9 cos — sm c; sni oq cos " 
sin 9 cos ^0 + cos 9 sin Oq cos 7 
^1 - (cos 9 cos 6*0 - sin 8 sin 6*0 cos 7)2 d^} 
sin 9 (sin 9 cos 9q + cos 9 sin 9q cos 7) 
1 — (cos 9 cos 6*0 — sin 9 sin 6*0 cos 7)2 ' 



^^^^ ^ l^uiii w ly^ 

1 — (cos 9 COS 6*0 
0(p sin 00 sin 7 i99 

9?/; 1 ~ (cos 9 cos 9q — sin 9 sin 6*0 cos 7)^ dtp ' 



(89) 

(90) 
(91) 



Notice that, by Eq. (|8^) and Eq. (|90D, the vector is proportional to sin 9 for all generic values of 0Oi except for 
(?o = OjTT. This means that, generically, the vector M2 vanishes at 9 = 0, which, by Eqs. (|8l|)-(]8^), represents source 
points {9, (j)) along the optical axis. These are the caustics of the past lightcone of the observer. With the metric 
Eq. ([Tq) we have 



Pil - 2V2ml) 




M2-M2 = - 

Ml -Ma = 



sin^9 



(92) 

(93) 
(94) 



Notice that these two connecting vectors are orthogonal along the lightcone. This is a feature of the spherical symmetry 
of the spacetime, which induces the axial symmetry of the lightcone. No two other Jacobi fields in general will be 
orthogonal, so these two are a special basis of Jacobi fields. However, they don't have unit length. In accordance with 
the previous subsection, we define two connecting vectors now that are orthonormal at the source's location: 



'-*(Mi-Mi) •\/-*(M2-M2) 
In our example, we will first examine a circular source, which we can express by 

Z"(/,t) = costM^{l)+smtM^{l). 



(95) 



(96) 



Because Mf and remain orthogonal, one can see that the circle at is distorted to an ellipse at any other point 
L with semiaxes 



and the ratio of the semiaxes of the ellipse at the observer's location is given simply by 



M2 


M2 


*(M2 


•M2) 


Ml 


Ml 


*(Mi 


•Ml) 



(97) 
(98) 



n 



/0(M2-M2 



_ *(Mi-Mi)"(M2-M2) 

o(Mi-Mi) " \} *(M2-M2)0(Mi-Mi)' 



(99) 



where we use the notation '^{v-w) =v-w\ig. With Eqs. (p2)-(93) we obtain 



n 



sin9(?o) de/dip\i. 




2V2mlp) 




- 2V2mh) 


99/9^1/0 sinOiQ \ 




2V2mlp) 




- 2V2mh) 



(100) 
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The factor gQ/g"^' \ ^ represents the ratio of the norms of the two starting connecting vectors at the observer's location 
(the apex of the lightcone) , at which all connecting vectors vanish. However the ratio of the norms is finite and can 
be calculated in closed form via a limiting procedure Q, obtaining 



The other factors in the ratio Eq. (100) must be evaluated numerically. Figure || shows a plot of the inverse ratio TZ 



at fixed value of the inverse radial distance to the source /*, as the image angle varies. We can see that the inverse 
ratio, representing the minor axis over the major axis, approaches 1 for large image angles (lightrays that pass far 
from the lens), but is everywhere smaller than 1, which means that the images are elongated along the polar direction 
around the optical axis in the lens plane (see Fig. The inverse ratio vanishes at the caustics, namely, when the 
source lies along the optical axis, behind the lens. This means that the distortion increases dramatically when the 
source is close to a caustic, leading to the formation of arcs. We have chosen to plot the inverse ratio rather than 
TZ itself, because the major axis becomes infinite for sources lying at the caustics, which makes TZ less convenient to 
plot. 

In the general case of the images of elliptical sources, they are given by the construction Eq. (p4). In our case 



(102) 



which shows, in the limit as I —^ Iq, that the magnification matrix, and therefore, the Jacobian matrix as well, in this 
case is symmetric, as is in the case of standard gravitational lensing. The stretching of the elliptical source occurs 
then along the eigendirections of J, which are the directions of Mf and AI2, since the matrix is diagonal. This means 
that if the source's semiaxes are aligned with the vectors (Mf , Mf ) then the image will not be rotated with respect to 
the source, so the distortion parameter measuring orientation distortion will vanish, Ds — 0. However, if the source's 
semiaxes are not aligned with (Af", Afj ), then the image will be rotated, and Ds will not vanish. See Fig. |^ 



IV. NON-PERTURBATIVE USE OF FERMAT'S PRINCIPLE IN LENSING 



In this section, we want to give a view of distortion that is based on the existence of families of null surfaces - 
this view in turn leads to a non-perturbative version of Fermat's principle applied to lensing. Although to apply 
it in practice is difficult, there being so few cases where enough null surfaces are known, it can be implemented 
perturbatively. In any case, it yields a slightly different picture of what is taking place in lensing. 

We begin by assuming that a two-point function, G{x°-,XqXX)i is known where for each value of the parameters 
(C,C): G = 0, represents a null surface containing the two points {x"',Xq). One should think of as being fixed 
and G = as representing a two-parameter family of null surfaces containing a;g; the parameters could be thought 
of as representing the sphere of null directions at Xq. The gradient of G, i.e., daG = Ga = £a, is a null covector, 
g'^'^GaGb = 0. 

Though we make no use of it now, we mention, for later use, that there is a special choice of a two point function 
that is often quite useful; if we begin with an asymptotically flat spacetime with the past null boundary, 5*^ x R, 
referred to as 3^, and coordinatized by (C,C) on the and u on the R, then the future lightcone of each point 
(m, C) of 3^ can be written as m = Z{x"', Ci 0- Since by definition the gradient of Z, i.e, daZ is a null covector, then 
by defining 

G(x-, xg, C, C) = Zix'\<:,0 - Z{x-oX,C), (103) 

we have the two point function satisfying our null surface conditions. 

We want to construct the envelope of G{x°',Xq,CiC) — imder variations of (C,C)- The envelope is a special 
null surface, namely the lightcone through Xq. Explicitly, this is done by setting to zero the (C,C) derivatives of 
G{x'^,Xq, C)- The envelope is then given by: 

G(x",a:S,C,C) = 0, 
dcG{x-,x-oX,C) = 0, 

5^G(x^a:S,C,C) = 0. (104) 

If the second two equations could be solved for (C, C) — (r(a;°, xg), r(a::", xg), they could be substituted into 
G{x^,Xq,(X) — Oj yielding g{x°',XQ) = G(x", a;g, F, F) = 0, the analytic expression for the lightcone of Xq. Whether 
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or not they could be solved, the three equations, Eqs. (104), could be solved for some three of the four x° in terms 
of the fourth coordinate x* and the parameters, i.e.. 



x'=a:Xx*,xS,C,C). (105) 



Equation (105) is the parametric form of the lightcone of xg, where the parameters, (C, C); label the generators of the 
lightcone (((^, C,) and play the same role as {9, 4>) did in previous sections.) 

It is from this parametric form of the lightcone that we will construct the Jacobian matrix of the lens equation. 
Note that if the x* = r \s a. radial type of coordinate, then x°(r, C) — {x* — a;*(r, .Tq, Cj C)j ^* = '''} form the 
lens and time of arrival equations where the observer is at Xq{t). The derivatives of with respect to (j', CiC) ^'^^ 
respectively the null geodesic tangent vector and the geodesic deviation vectors. 

We first point out that the construction we just gave can be given an alternative meaning; if the equation 
G{x°', XqXX) =0 is solved for the coordinate i = a;°, so that 

G = ^ t = G(2;",2;S,,0, 

then the envelope construction is identical to "extremizing" the time, t, with respect to (^, Q) variations. We thus have 
constructed the lightcone of via Fermat's principle. It is then clear that the construction of the lens equation that 
was given in Sec. II via the Fermat potential is a special case of the more general construction we have just outlined. 
It will be seen that from this construction there are natural coordinates so that the Jacobian matrix is symmetric. 



By treating the three equations, Eqs. (104) and x* = r = constant as four implicit equations for the determination 



of as functions of (C, C) we can, by differentiation with respect to (C, C) construct the deviation vectors, M°- — d^x"" 



and M = d^x"^ . Explicitly, the derivative of Eqs. (104) and x* = r with respect to C, are 



Arc, + acG(a;^ xf;, C, C) = M'^G, = 0, 
M'^acGa(x^xg,C,C) + ac9cG(a^^xg,C,C)-0, 

M-d^Ga (x^ , c, + ac9^G(a:^ , C, C) = 0, 

M^Sl = [dcx*) = 0, (106) 

and the complex conjugate equations are 

M°G„ + a^G(x^ C, C) = ^"Ga = 0, 
M"acGa(a;", xg, C, C) + c^9cG(x^ xg, C, C) = 0, 

M\Ga(x^ xg, CO + a^c^G(x^ xg, c, c) = o, 

= 0. (107) 

In the first of each of the two sets we simplified using the envelope condition. We thus see that there are only two 
non-vanishing tetrad components for both Af ° and M . We define these components as 

M+ = M"9cGa(x",a;S,C,C) - -d^d(^G{x- .x'^.XX), 
M- = M-d^a{x\xl,U)^-&^dcG{x%xlU). 
It = M''dcGa{x%x-oXX) = -c%9cG(x^xS,C,C), 
= M\Ga (x^ , C, C) = -d^d^Gix- , , C, , 



and we have the tetrad version ^ of the Jacobian matrix, J, i.e., 

(108) 




-d^dcG -a^d^G 
-d^d^G -d^d^G 

(Note that d^Ga and d^Ga play the same role as did in previous sections). We see that the Jacobian matrix 
is symmetric when expressed in the tetrad basis. A local coordinate system (depending on the value of (C, C) or 
equivalently on the relevant null geodesic between source and observer) can always be introduced so that the tetrad 
components of 3 are also the coordinate components. 
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It should be noted that though there is no easy way to relate Z directly to the shear and divergence of the 
lightcone congruence, there nevertheless is a rather startling relationship between ^ and certain asymptotic quantities 
(asymptotic s hear s) in the case of asymptotically flat spacetimes when the two-point null surface functions G are 
given by Eq. ( |103| ). This relationship is both startling and (at least so- far) rather incomprehensible. At the moment 
it does not se em t o have any bearing on astrophysically measurable quantitites and is only of theoretical interest. 



Using Eq. (|103D , Eq. ( |108|) becomes 

a(x^xg,c,c) 



- [d^dc z-d^dcZo] - [d^d^z - d^d^Zo] 



Though, as we mentioned earlier, the function u = Z{x°' , represents the future lightcone from the point (u, C, 
of 3^ , it has an alternative reciprocal meaning: If is held constant and (C,C) ^re varied over the sphere, then 
u ~ Z{x°' , C,, Q describes the the intersection of the past lightcone of the point x"" with 3~, the past lightcone cut of 
the point x"" . This past lightcone, as it approaches in the limit 3^, has an asymptotic shear (Jx{x°' tCX) 

It has been shown in the completely different context of the Null Surface Reformulation of general relativity [p3|j24[ 
that ^(a;'', Xg, C) has an asymptotic interpretation, namely 



Z[x'',xlXX) 



-[a,(a:^C,C)-^.o(^g,C,C)] -[R-Ro] 

-[R-Ro] - px (x^ c, c) - i^tx.o] 



where R and Rq are, respectively, measures of the curvature of each of the lightcone cuts from a;° and Xq. This matrix 
constructed from the differences between the two asymptotic shears are surprisingly a measure of the changes in 
geodesic deviation of geodesies going from a;° and Xq. It is this fact that we find both startling and incomprehensible. 



V. CONCLUSIONS 



We have developed a non-perturbative description of distortion of images to be used in a generic spacetime in a 
manner analogous to the case of standard gravitational lensing, where a set of simplifying assumptions is normally 
used. 

We have introduced three shape parameters which describe a small, elliptically shaped cross-section of a pencil of 
rays from an "elliptical" source in the past lightcone of an observer. These are the cross-sectional area A of the pencil 
(or its angular analog fi); the ratio of the major to minor semiaxes of the ellipse, TZ] and the orientation of the ellipse, 
5 with respect to some arbitrary reference axis. By propagating these quantities towards the observer, we have found 
the shape parameters of the image, which lies on the observer's celestial sphere. Based on the comparison of the 
shape parameters of the source and image, we have defined three distortion parameters, the values of which provide 
a quantitative measure of the severity of the distortion introduced by the gravitational field on the path of the pencil 
of rays. 

Our work in this paper is based on the use of Jacobi fields to propagate the source along the lightcone. Our 
approach can be directly related to astrophysical gravitational lensing in the weak-field regime where the thin lens 
approximation is used by noticing that, within such a regime, the Jacobian matrix A is equal to our Jacobian matrix 
*J scaled by the distance to the source, s*, namely: 




When the weak-field thin-lens regime is not imposed, however, there is no geometric meaning to J / s. This is a reason 
why we prefer to refer to J itself as the Jacobian and consider it as the closest meaningful analog of the matrix A in 
use in astrophysical gravitational lensing. 

Although our calculations in this paper do reproduce the formalism of image distortion of standard lensing theory 
in the regimes of weak fields, thin lenses and small angles, they do not make any assumptions about the strength 
of the gravitational field or about the thickness of the mass distribution. Interesting prospective applications of 
this approach include the development of alternative approximation techniques to the lensing problem in which the 
gravitational field is weak (so that the metric can be found from the linearized Einstein equations, or equivalently, by 
adding Schwarzschild-like contributions) but where one utilizes a fully three-dimensional mass distribution instead of 
multiple lens planes. 

Having defined the shape parameters of the pencil of rays in terms of the Jacobi fields of the observer's lightcone, 
our next goal is to clarify the relationship between the optical scalars of the lightcone and the change in the shape 
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parameters of the pencil of rays. We derive this relationship in our companion paper [gj. In addition, in our companion 
paper |^ we derive the expression of the Jacobian matrix A in terms of the optical scalars of the lightcone in the 
weak-field thin-lens approximation. 
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o 




FIG. 1. The lensing problem in standard gravitational lensing. The point O represents the observer. The point S represents 
the source. The lens mapping can be interpreted as a map that takes the point / in the lens plane into the point S in the 
source plane. 
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observer' s\ 
worldline 



source' s 
worldtube 




FIG. 2. Diagram of image distortion. The source's worltube intersects the observer's past lightcone in a region free of 
caustics. The source's visible shape is defined by the intersection. The pencil of rays between the source's shape and the 

observer carries the shape of the source into the shape of the image, on the observer's celestial sphere. The pencil of rays can 
be described by geodesic deviation vectors (connecting vectors of the observer's lightcone) from a central null ray connecting 
the center of the source's shape to the observer. 
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FIG. 3. Vectors associated with the null geodesic that joins the source to the observer. The vector I is tangent to the null 
geodesic. The vectors (61,62) are spacelike, orthogonal to ^, orthonormal and parallel progagated along the geodesies. The 
vectors (Mi,M2) are linearly independent Jacobi fields which coincide with (61,62) at the location of the source (not shown). 
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ce 





FIG. 4. The connecting vectors in Schwarzschild spacetime. The vector Mi = d/dij) Ues on the plane defined by the source, 
the observer and the lens. The vector M2 = d/d^f lies on the direction orthogonal to such a plane. The inverse radial distance 
to the source {I = (\/2r)~^) is used as a parameter along the null geodesic that connects the source to the observer. The 
function Q(l,ip) represents the angular location of the source with respect to the lens. 
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FIG. 5. Distortion by a Schwarzschild black hole. The inverse ratio of the semiaxes of the image divided by the inverse ratio 
of the semiaxes of the source is plotted for fixed distance between the source and the deflector, equal to the distance between 
the lens and the observer, as a function of the image angle. For large image angle the distortion decreases, since the ratio 
tends to 1 as V ^ -k 12. For small image angle, the geodesies that arrive at the observer go through a strong gravitational field, 
and the distortion is large. The distortion is infinite (where the inverse ratio is zero) for the images of sources that lie on the 
caustics the observer's lightcone. 
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Source plane 

FIG. 6. Distortion by a Schwarzschild black hole. The axial symmetry in the source plane leads to no orientation distortion 
if the semiaxes of the elliptical source are aligned with the connecting vectors Mf and . points in the polar angular 
direction, whereas M2 points along the radial direction on the source plane. Thus, only sources aligned with the radial direction 
experience no distortion of orientation. 
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